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^ ■ 1 Introduction 

^ : 

I The further progress in understanding the phase structure of Quantum Chromodynamics (QCD) 
at non-zero temperature and density is very important in view of the recent experiments at RHIC 
and also planned at LHC and FAIR. From the theoretical point of view, it can be divided into 
some regions of "phase space": perturbative QCD, lattice QCD, phenomenological models, etc. 

As known, perturbative QCD works well at asymptotically high densities [Tj or temperatures 
[2], where the QCD coupling ag is small. For lower densities and temperatures the MIT bag model 
can be of use [3]. It is a phenomenological model that effectively includes strong interactions via 
I a bag constant. 

The bag model has taken an interesting turn with the detailed analysis [3] of the lattice data 

[5]. The surprise of this analysis is that at zero chemical potential it reveals a term quadratic in 

temperature as the leading correction to the ideal gas term in the pressure. The most recent data 

[6] indicate that this is generic: that the equation of state of the MIT bag model gets modified 
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p{T) = aT^ - T^T^ - B , with T^^^ < T < T_ . (1) 

As noted, the novelty is the term, while the remaining two are the standard bag terms, with 
B a bag constant and a a parameter. A common choice is to take a from perturbation theory 
up to one loop order. r,„i„ is close to a critical temperature Tc (or some approximate "Tc" for 
a crossover). A small difference between Tc and T^i^ may vary with the model. T^^^ is set by 

*Also at Landau Institute for Theoretical Physics, Moscow. 
^Rob Pisarski called it a "fuzzy" bag model for the pressure [1]. 
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perturbation theory such that to leading orders it is applicable only for temperatures higher than 
T 

max • 

While numerical simulations on the lattice can be of use at non-zero temperature when the 
quark density is quite small, standard Monte Carlo techniques are not of use in cold dense matter 
because of poor convergence called the sign problem. On the other hand, there are arguments 
in the literature [7] that in QCD long perturbative series (or the UV renormalons) result in the 
so-called quadratic corrections. From this point of view, the term in ([T]) is nothing else but an 
example of the quadratic correction. If so, then it is natural to expect that at zero temperature 
the equation of state of the MIT bag model can also get modified by a quadratic correction a^ 

= - ^ilp^ - B , with < /i < , (2) 

where again 6 is a parameter to be fixed from perturbation theory up to one loop order. The 
/i^ term is a quadratic correction. /Xn^i^ is expected to be close to a critical value of while 
/iniax is set by leading orders of perturbation theory. Note that /i stands for the baryon chemical 
potential, here and below. 

Until recently, the lattice formulation and effective field theories were the main computational 
tools to deal with nonweakly coupled gauge theories. The situation changed drastically with 
Maldacena duality (AdS/CFT) [Tl] that resumed interest in another tool, string theory. The 
original duality was for conformal theories, but various perturbations (deformations) produce 
gauge/string duals with a mass gap, confinement, and chiral symmetry breaking |12j . 

In this Letter we continue a series of recent studies [131 El ES] devoted to a search for 
an effective string theory description of strong interactions. Since precise recipes for finding the 
string theory dual to QCD are still unknown, our strategy is based on deformations of AdS/CFT. 
The deformation we are pursuing turned out to be successful in providing a systematic approach 
to the quadratic corrections. Indeed, in [13], the quadratic correction was found in the two- 
current correlator. Later, the model was extended to Euclidean signature for computing the 
heavy quark potential, where the quadratic correction occurs as a liner term in the potential 
at short distances |14]. Subsequent comparison |16j with the meson spectrum made it clear 
that the model should be taken seriously. Moreover, it was also extended to finite temperature. 
As a result, the term in the pressure ([TJ was found [T5j. In addition, this model results in 
the spatial string tension [T^ and the expectation value of the Polyakov loop [18] which are 
remarkably consistent with the lattice. Thus, there are reasons to believe that this deformation 
also provides a good approximation for a string dual to cold quark matter. 



2 The Model 

Let us first explain the model to be considered. We take the following ansatz for the 10- 
dimensional background geometry which is a deformed product of the extremal Reissner- Nordstrom 
black hole in Euclidean AdSs and a 5-dimensional sphere (compact space xjl 



^Such a parametrization of the quark matter equation of state was first proposed in tiie context of phenomenol- 
ogy of neutron starts with quark interiors [8^ and further discussed in [91 110). 

''For completeness, we include a brief summary of the relevant results concerning the Reissner-Nordstrom black 
holes in the Appendix. 
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H = e(V2)-^ , / = (1 _ [z/z^ff{l + 2{z/z^f) , 

where z+ = {2/q^)&. The deformation is due to the same z-dependent factor H as those of 
|13^ [MKTS], with c being a parameter whose value will be fixed shortly. Note that ([3]) is smooth 
and complete at z = z+ such that the inverse period of t is equal to /3~^ = T = 0, with T the 
temperature. We also take a constant dilaton and discard other (if any) background fields. 

Given the background metric ([3]), we can now find the corresponding gauge potential from 
the condition of Weyl invariance on a string world-sheet. To leading order in a' it is given by 

= V^F,,, + O(a') = 0. (4) 

Here is in fact a renormalization group beta function on the world-sheet. 
For a pure electric potential Aq{z), becomes 

5.(^/55°V^5,^o) =0. (5) 

The solution is given by Aq{z) = Cie'^^^ + C2, with Cj constants. If we choose the constants so 
that at c = the solution is reduced to that of Reissner-Nordstrom (|A.6p . we find 




Finally, we impose the condition Aq{z^) = and as a result get // as a function of q 

M,) = ff(e"i-l). (7) 

Following AdS / CFT dictionary , we identify the parameters fj, and p, as defined in ()A.7p , 
with the baryon chemical potential and the baryon number density, respectively. 



3 Estimate of the Quadratic Correction 

Our first goal will be to analyze the non-deformed model. That is, we take c = for ([3]) and ([6|). 
Using the formulas and (|A.7p . we can show that the equation ([7]) yields, up to a constant 
multiple, a unique solution p{fi). Explicitly, it is given by 

pip) = 4bp^. (8) 

This is the desired result, corresponding in QCD to the fact that for low temperatures and 
large chemical baryon potentials the baryon number density is proportional to the cube of the 
potential. 

To fix the constant of proportionality, we need some knowledge of the exact string theory dual 
to QCD or some additional insight. Since the former is beyond our grasp at present, we match 
the parameter b with that of perturbative QCD neglecting perturbative interactions among the 
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quarks. In doing so, we first find tlie pressure by integrating ^ = P- In terms of tlie quark 
chemical potential fig = ii/Nc, it is p{pq) = bN^fig. Finally, we have 

(9) 



127r2 m ' 



Here Nf is the number of quark flavors and Nc is the number of colors. 

Now let us return and discuss the deformed model. At large baryon density (or equivalently 
at large q) it is reasonable to represent ([7]) as a series Yln=i T^i^)^ ■ take the two 

leading terms of the series, then we can easily invert the function p{q). Finally, using the equation 
()A.7p . we find the leading correction to ([8]) 

p(/x) = 46(/i3 -1^2^ + 0(1)) , with ^^* = ^]|l■ (10) 

In the homogeneous case the pressure is obtained by integrating the above expression over fj, 

p(^) = 6(/-/.V + 0(l)). (11) 

This is our main result. It includes the /i^ term, as expected. 

Making an estimate requires some numerics. First, let us consider the light {u, d) quarks. In 
this case, the value of c is fixed from the slope of the Regge trajectory of p{n) mesons [19J. This 
gives So, we find for the value of the quadratic correction 

pl^A.lGfN'^. (12) 

In contrast, a simple estimate of the corresponding coefficient of perturbative QCD with Nf = 2 
results in [1] 

^(m2+m2)«6xlO-4GeV2. (13) 

Here we have used that iriu = 3 MeV and rud = 6 MeV. 

Thus, our model predicts that the p"^ term being negligible in a pure perturbative region 
A*max < P gets strongly enhanced in the intermediate region p^^^ < p < /i^ax- 

Next, let us discuss the effect of the strange quark. For Nj = 3, (fT3|) becomes 9(m2 +m^+m2) 
which is certainly valid near the upper limit p^n^x, where p 3> 3ms- A simple algebra shows that 
its value is of order 0.1 GeV^, with ~ 0.1 GeV. It is still smaller than p2|) . so the effect of 
the strange quark is not dominant. 

Finally, let us estimate the range of fi for the model of interest. A crude estimate of the 
lower limit can be made by using the positivity of the baryon density and the pressure. It gives 
that //min is of order If we assume that like at finite T on the lattice [6], where T^i^ ~ l-5Tc, 
in the model of interest fi^i^ ~ 1.5/Uci then using ([12]) we arrive at a critical chemical potential 
of 1.3 GeV which is reasonable phenomenologically. A crude estimate of the upper limit can be 
made by assuming that at fi = /^max the contribution of the IJ? term in the pressure is one order 
of magnitude smaller than that of the leading //^ term. This gives /U^ax ~ 3.3/U^ or, in terms of 
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4 Concluding Comments 



(i) Having derived the equation of state, we can easily develop finite /i thermodynamics. In 
particular, for the energy density, we have e = 6(3/x^ — /x^/i^ + 0(1)). Combining with (fTTj) . we 
find the expression for the trace anomaly 

i^ = 264 + 0(l)- (14) 
In addition, for the speed of sound C| = ^ we get 

C2(/.) = i(l-l|| + 0(l)). (15) 

All the above formulas are similar to those of fT5] at finite T. 

(ii) Here we used the model based on the deformation of the Reissner-Nordstrom solution. 
Certainly, such a phenomenologically motivated way is out of the mainstream of (academic) 
AdS/CFT, where the background geometry follows from the equations of supergravity and funda- 
mental matter is introduced via D-brane embeddings in the probe approximation with Nc ^ A'^^. 
One of the advantages of our approach is that it allows us to incorporate the backreaction due 
to the gauge potential on the background geometry. What really fits better to QCD remains to 
be seen. 

(iii) In the phenomenological parametrization of [9] the coefficient in front of the /i^ term arises 
from the strange quark mass as well as color superconductivity. As a result, it is proportional to 

— 4A^. Its value is one order of magnitude smaller than ours (jl2p . 

The formula ([2]) was also suggested, by analogy with the deformed bag model in the 
context of Quarkyonic matter [10]. Our interpretation of the ^'^ term as a power correction 
differs from that of [10], where it is interpreted as due to nonperturbative corrections. However, 
in the intermediate region of interest some matching conditions between the two regimes may be 
possible. 
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Appendix 

In this appendix we review the relevant results concerning the Reissner-Nordstrom solutions in 
five dimensions. Most of this material can be found in [20]. 

For the Einstein-Maxwell action with cosmological constant, we take 

Here Gn is the 5-dimensional Newton's constant. 
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With a pure electric gauge potential 



^o(r)=^(-^^+/x , (A.2) 



21 r2 

a solution of the equations of motion for the metric (with Euclidean signature) takes the spheri- 
cally symmetric form 

2 2 

ds^ = fde + r^dr^ + r^d^l , / = l- !^ + ^ + ^. (A.3) 
The parameters m and q are respectively related to the mass and charge of the black hole as 

3Vol(S3) ^ \/3Vol(S3) 

Here Vol(S^) is the volume of a unit 3-sphere. 

The solution (jA.Sp is asymptotic at r = oo to x S^. A scaling that reduces it to a solution 
with X may be made as follows. If we introduce a dimensionless parameter A and make 
the transformation r — )• X^r ,t — )• X~it,m — >• Xl^m,q — )■ X^l^q, then in the large A limit we 
obtain 

/2 

ds'^ = — [fdt^ + f-^dz^ + dx^) , / = 1 - mz^ + , (A.5) 

where z = P/r. In the process, we have also introduced local coordinates yi near a point P G S'^ 
such that dQ,'^ = ^ dy^^, and then set Xi = X^lyi. 

Having derived the desired solution for the metric, we can easily obtain that for the gauge 
potential. From (|A.2p . we have 

Aoiz) = i(^-^qz^ + f?j. (A.6) 

When we go to R'^ x S^, we get that the radius of is proportional to A 4 and so diverges 
for A — )• oo. Hence, the corresponding volume is also becoming infinite and looks like V3 = 
A4^^Vol(S^). If we introduce the charge density p = Q/V3, then the second equation of ()A.4p 
becomes 

p = sVSbq , (A.7) 

where b = / {IIttGn). The difference between S'^ x S"^ and R^ x S"*^ is obvious: in the first case 
q is related to the charge of the black hole, while in the second case it is related to its charge 
density. 

The metric (jA.Sh is smooth and complete if the period of t is /3 = tj^£^, where z+ is the 
smallest real positive root of f{z) = 0. 

For T = 0, the black hole becomes extremal so that 4m^ = 27(/^. In this case the function 
f{z) takes the form 

f=[l-{z/z+f)\l + 2{z/z+f), with z+ = (2/g2)i. (A.8) 



The parameter ^ is reserved for future use. 
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